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I.  INTRODUCTION 


Progress  has  been  achieved  during  1989  and  1990  in  the  study  of 
nonlinear  systems  which  generate  chaotic  1/f  fluctuations,  in  the 
application  of  the  Quantum  1/f  Theory1’3  to  various  materials  used  in 
small  and  ultrasmall  electronic  devices,  and  in  the  application  of  the 
Quantum  1/f  Theory  to  electronic  devices. 

The  close  similarity  cf  the  classical  and  quantum  1/f  tneories,  ana 
the  initial  development  of  the  quantum  1/f  theory  by  the  author  out  of  his 
efforts  to  quantize  his  classical  turbulence  theory,  have  led  to  sustained 
efforts  of  the  author  aimed  to  integrate  all  his  theories  as  various  forms 
or  realizations  of  a  fundamental  notion  of  chaos  in  nonlinear  systems. 
During  this  grant  period,  these  efforts  finally  beared  fruits.  A  general 
sufficient  criterion  was  formulated,  allowing  to  identify  the  nonlinear 
systems  which  exhibit  1/f  spectra.  This  criterion  is  presented  for  the 
first  time  in  Sec.  II  below.  It  is  followed  in  Sec.  Ill  by  examples,  in  which 
the  criterion  is  applied  to  the  classical  and  quantum  mechanical  forms  of 
the  author's  1/f  noise  theory.  These  examples  clarify  the  physical 
meaning  of  the  new  criterion. 

For  the  practical  application  of  the  Quantum  1/f  Theory  it  is 
necessary  to  derive  the  quantum  1/f  fluctuations  of  various  kinetic 
(transport)  coefficients  which  characterize  the  materials  used  in 
electronic  and  microelectronic  applications,  from  the  author's 
fundamental  quantum  1/f  formula.  The  latter  is  applicable  only  to  cross 
sections  and  rates  of  elementary  processes.  Most  important  is  the 
calculation  of  mobility  fluctuations  in  Si,  GaAs  and  Hgi-xCdxTe.  An 
earlier  calculation  (Kousik,  Van  Vliet,  Handel,  1985)  of  mobility 
fluctuations  in  Si  and  GaAs  is  replaced  in  Sec.  IV  by  a  more  rigorous 
calculation,  based  on  the  new  quantum  1/f  cross-correlations,  developed 
under  the  previous  AFOSR  Grant,  and  presented  in  the  Final  Technical 
Report  AFOSR  -85-0130.  The  new  calculation  yields  increased  1/f  noise, 
and  is  in  vey  good  agreement  with  the  experiment.  Although  not  presented 
here  yet,  we  also  performed  a  Monte  Carlo  simulation  for  HgCdTe  during 
this  period.  The  simulation  has  yet  to  be  improved  and  compared  with  the 
experiment. 
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Finally,  we  have  tried  to  improve  the  application  of  quantum  1/f 
theory  to  the  collector  noise  of  bipolar  transistors.  This  short  calculation 
is  presented  in  Sec.  V. 

In  the  same  time,  many  new  contributions  to  the  quantum  1/f  theory 
and  experiment  were  published  by  other  workers  in  the  field,  considerably 
advancing  the  field  of  infra-quantum  physics  and  quantum  1/f  noise  in 
high-technology  applications.  These  new  contributions,  as  well  as  new 
PhD  thesis  work  in  this  field  and  contributions  presented  at  the  IV 
Conference  on  Quantum  1/f  Noise  and  other  Low-Frequency  Fluctuations, 
are  included  in  the  updated  General  Quantum  1/f  Bibliography  appended  to 
this  Report. 


II.  A  SUFFICIENT  CRITERION  FOR  1/F  NOISE  IN  NONLINEAR  SYSTEMS 

Consider  a  n-dimensional  nonlinear  system  described  in  terms  of  the 
dimensionless  function  Y(x,t)  by  the  mth  order  nonlinear  dynamical 
equation 

dY/dt  +  F(x,  Y,  dY/dx-)  ...dY/dxn,  d2Y/dx-|2 . d™Y/dxnm)  =  0  (1 ) 

I  f 

F[Xx,  Y,  dY/(Xdx1)...dY/(^dxn),  d2Y/(Mx-|)2 . dmY/(^dxn)m] 

=  XPF(x,  Y,  dY/dxi„. dY/dxn,  d2Y/dx-|2 dmY/dxnm)  (2) 

for  any  real  number  X,  Eq.  (1)  is  said  to  be  homogeneous.  Performing  a 
Fourier  transformation  with  respect  to  the  vector  x,  we  get  in  terms  of 
the  Fourier-transformed  wavevector  k  the  nonlinear  integro-differential 
equation 

dy(k,t)/dt  +  G[k,  y(k,t),  k-|  y(k,t)...kny(k,t),  k-|2y(k,t) . knmy(k,t)]  =  0,  (3) 

where  y(k,t)  is  the  Fourier  transform  of  Y(x,t).  Due  to  Eq.  (2),  G  satisfies 
the  relation 
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G[Xk,  y,  Xk^y.-AknY,  (Xki)2y . (kkn)my] 

=  XPG[k,  y,  k-j  y...kny,  k-(  2y . knmy],  (4) 


Eq.  (3)  can  thus  be  rewritten  in  the  form 

dy/d(t/A.P)  +  G[Xk,  y,  M-|y...Xkny,  (Xk-j )2y . (Xkn)my]  =  0,  (5) 

Taking  A.=1/k,  where  k=|k|=(x-j  2+....+xn2)",//2,  and  setting  kPt=z,  we  notice 
that  k  has  been  eliminated  from  the  dynamical  equation,  and  only  k/k  is 
left.  This  means  that  there  is  no  privileged  scale  left  for  the  system  in  x 
or  k  space,  other  than  the  scale  defined  by  the  given  time  t,  and  expressed 
by  the  dependence  on  z.  We  call  this  property  of  the  dynamical  system 
"sliding-scale  invariance 

In  certain  conditions,  instabilities  of  a  solution  of  Eq.  (1)  may  generate 
chaos,  or  turbulence.  In  a  sufficiently  large  system  described  by  the  local 
dynamical  equation  (1),  in  which  the  boundary  conditions  become 
immaterial,  homogeneous,  isotropic  turbulence,  (chaos)  can  be  obtained, 
with  a  spectral  density  determined  only  by  Eq.  (1).  The  stationary 
autocorrelation  function  A(x)  is  defined  as  an  average  over  the  turbulent 
ensemble 

A(x)  =  <Y(x,t)Y(x,t+x)>  =  J<y(k,t)y(k,t+x)>dnk  =  Ju(k,z)dnk  (6) 

Here  we  have  introduced  the  scalar 

u(k,z)  =  <y(k,t)y(k,t+x)>  (7) 

of  homogeneous,  isotropic  chaos  (turbulence),  which  depends  only  on  |k| 
and  z=kPx.  All  integrals  are  from  minus  infinity  to  plus  infinity.  The 
chain  of  integro-differential  equations  for  the  correlation  functions  of 
any  order  obeys  the  same  sliding-scale  invariance  which  we  have  noticed 
in  the  fundamental  dynamical  equation  above.  Therefore,  in  isotropic, 
homogeneous,  conditions,  u  can  only  depend  on  k  and  z.  Furthermore,  the 
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direct  dependence  on  k  must  reflect  this  sliding-scale  invariance,  and  is 
therefore  of  the  form 

u(k,z)  =  k_nv(z).  (8) 

Indeed,  only  this  form  insures  that  u(k,z)dnk  and  therefore  also  the 
corresponding  integrals  and  multiple  convolutions  in  k  space  have  the 
necessary  sliding-scale  invariance. 

According  to  the  Wiener-Khintchine  theorem,  the  spectral  density  is 
the  Fourier-transform  of  A{t), 

Sy (f)  =  Je2rcifTA(T)dT  =  (1/f)  Je2rcit,Jk,-nv(z)dnk'dt'  =  C/f,  (9) 

where  we  have  set  fc=t',  kn=fk'n,  z=kn<t=k'ntl,  and  the  integral 

0  =  Je2rcit,Jk'-nv(z)dnk'cit'  =  Je2^it,Jk"-nv(k,,n)dnk,,dt,  (10) 

is  independent  of  f.  We  have  defined  the  vector  kM=t'1/n  k. 

In  conclusion,  we  have  shown  that  if  the  equation 

dY/dt  +  F(x,  Y,  dY/dx-|  ...dY/dxn,  d2Y/dx-| 2 . dmY/dxnm)  =  0  (11) 

with 

F[Xx,  Y,  dY/(X.dx1)...dY/(Mx)n,  d2Y/(Xdx-j)2 . d™Y/(Xdxn)m]  =  XPF(x,  Y, 

dY/dx-j  ...dY/dxn,  d2Y/dx12 d^Y/dx^)  (12) 

admits,  in  the  limit  of  weak  dissipation,  quasistationary  homogeneous 
isotropic  chaotic  (turbulent)  solutions  which  are  practically  independent 
of  the  nature  of  the  instabilities  or  bifurcations  (or  even  stirring  forces) 
which  have  caused  the  chaotic  state,  the  corresponding  spectral  density 
must  be  proportional  to  1/f. 

We  note  that  the  solution  (7)  leads  to  a  weak  (logarithmic)  divergence 
of  the  integral  over  k  in  the  last  form  of  Eq.  (6)  and  in  Eq.  (9).  This  seems 
to  contradict  the  fact  that  in  practical  applications  the  autocorrelation 
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function  A(t)  is  finite,  and  its  value  at  t=0  is  usually  given  in  the  problem 
at  hand  However,  in  practice  one  never  deals  with  an  infinite  volume,  and 
the  physical  wave-vectors  are  also  limited.  For  instance  in  fluid 
dynamics,  wave  vectors  exceeding  the  reciprocal  average  distance 
between  neighboring  fluid  molecules  correspond  to  thermal  motions,  and 
are  therefore  no  longer  meaningful  for  the  hydrodynamic  treatment.  Due 
to  its  logarithmic  character,  the  divergnce  is  thus  without  practical 
importance.  Nevertheless,  for  a  given  level  of  chaos  A(0),  we  can 
construct  an  approximate  solution 

u(k,z)=kr_nv(z),  (13) 

with  0<r«1,  which  avoids  the  divergence  at  k=0.  To  get  the  correct  chaos 
level  with  k<k0,  k0  being  an  upper  cutoff,  we  set 

u  =  r[A(0)/v(0)]ko-rkr'nv(z).  (14) 

This  yields  for  t=0  the  result  A(0)  when  we  integrate  over  dnk  with  an 
upper  limit  k0.  We  notice  that  r  is  present  both  as  a  general  factor,  and  as 

a  small  defect  in  the  exponent  of  k.  This  is  a  general  feature,  present  both 
in  classical  and  quantum  nonlinear  systems  with  1/f  noise.  In  the  limit 
r->0,  the  approximate  solution  tends  to  become  exact.  In  the  classical 
homogeneous,  isotropic,  turbulence  theory,  r  can  be  arbitrarily  small, 
while  in  the  quantum  1/f  theory  (quantum  electrodynamics),  r  =  aA  = 
(2a/37t)(Av/c)2  «  1/137  is  the  well-defined  infrared  exponent  of  the 
process,  with  a=e2/hc=1/137,  as  we  know  from  the  theory  of  infrared 
radiative  corrections. 
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III.  EXAMPLES 


The  general  criterion  developed  in  the  preceeding  section  will  now 
be  illustrated  on  the  basis  of  some  examples. 

111.1  Classical  Turbulence  Theory  for  the  Current 
Carriers  in  Semiconductors 

In  the  case  of  homogeneous,  isotropic  turbulence4'6  caused  in  the 
electron-hole  plasma  of  an  infinite  sample  of  a  symmetric  intrinsic 
semiconductor  by  dynamical  instabilities  of  any  kind,  we  start  from  the 


equations 

vv+  =  (e/2c)vxB  -  (1/n)VP,  (15) 

w  =  2e[E  +  v+xB/c]  -  (2/n)V(Pp  -  Pn),  (16) 

V  v+  =  0  (n  =  const),  (17) 

V  xE  =  -(1/c)9B/3t,  (18) 

V  xB  =  27tenv-/c,  (19) 

V  B  =  0.  (20) 


Here  n  is  the  total  carrier  concentration  including  an  equal  number 
of  electrons  and  holes,  v/e  their  reciprocal  mobility  assumed  to  be  the 
same  for  electrons  and  holes,  Pn  and  Pp  the  partial  pressures  of  electrons 
and  holes,  P  the  total  carrier  pressure,  2v+'  v*  the  sum  and  the  difference 
of  the  carrier  drift  velocities.  Inertial  terms  proportional  to  the 
effective  masses  of  the  carriers,  as  well  as  electrostatic  terms  and 
compressibility  terms  have  been  neglected  here  in  a  consistent  way4'7, 
because  we  are  interested  in  the  low-frequency  domain  only.  Although  we 
do  not  work  this  out  here,  this  system  of  equations  can  be  shown  to  admit 
an  energy  theorem.  Performing  a  Fourier  expansion,  we  obtain 
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vv  +  (k)  =  (e/2c)Zk*v'(k’)xB (k-k’)  -  (i/n)k-P(k),  (21) 

vv-(k)  =2e[E(k)+  Ik'V  +  (k')xB(k-k,)/c]-(2i/n)k(Pp-Pn),  (22) 

kv  +  (k)  =  0.  (23) 

ik  xE(k)  =  -(1/c)3B(k)/at,  (24) 

ik  xB(k)  =  27tenv-(k)/c,  (25) 

k-B(k)  =  0,  (26) 

Substituting  E  from  Eq.  (22)  into  Eq.  (24),  we  obtain  with  Eq.  (25) 

3B(k)/3t  +  fik2B(k)  =  ikx  £«•  v+(k')  xB(k-k’),  (27) 

where  p  =  c2v/47ine2.  Eqs.  (21)  and  (25)  yield 

v+(k)  =  (i/47cvn)  Zk«{B(k")[k"  B(k-k")]  -  k"[B(k")  B(k-k”)] 

(28) 


-(k/k2)(1  -5ki0)[k"  B (k-k")  k-B(k")  -  B  (k")  B(k-k")  (k  k")]}  . 

Substituting  this  into  Eq.  (27),  we  obtain  the  fundamental  dynamical 
field-equation  of  turbulence  in  ihe  electron-hole  plasma  of  a  symmetric 
intrinsic  semiconductor 

dbp(k.t)/dt  +  fik2bp(k,t)  =  Zk'kH  bj(k-k',t)  b|(k",t)  bm(k’-kH,t) 

(29) 

•(kj5ps-ks5pj)[ks5|m-km"5|s  +  (ks7k’2)(1  -5k-io)(km,,k|’-k’k"5|m)]> 

in  terms  of  b  =  B/V27tvn.  This  dynamical  equation  has  the  form  of  Eq.  (3), 
with  p=2  in  Eq.  (4),  and  with  G  defined  as  the  r.h.s.  minus  the  term  in  k2  on 
the  l.h.s..  Our  sufficient  criterion  thus  tells  us  that  this  nonlinear  system 
will  yield  a  1/f  spectrum.  We  present  below  -he  direct  derivation  for  this 
example. 


Multiplying  Eq.  (29)  with  ba*(k,t-x)  and  taking  the  average  over  a 
statistical  ensemble  which  represents  our  notion  of  stationary 
turbulence,  we  obtain  in  quasi-stationary  conditions 

dwap(k,t)/3x  +  M-k2wap(k,x) 

=  Xk,k"<ba*(k>t-x)bj(k-k')t)b|(kM,t)bm(k,-k")t)>Rjmip.  (30) 

with  wap(k,x)  =  (Ly27i:)3<ba*(k,t-x)bp(k,t)>,  where  L  is  the  edge  of  the  cubic 
normalization  box,  and 


Rjmip*(kj6ps-ks5pj)[ks8|m-km'6ls  +(ks7k'2)(1  -5k'>o)(km"k|,-k'k"5|m)].  (31) 

Multiplying  Eq.  (29)  with  more  magnetic  field  components  and 
averaging,  we  obtain  equations  connecting  the  fourth-order  correlation 
tensor  to  the  sixth-order  tensor,  and  so  on4-7.  To  end  this  infinite  chain 
of  equations  for  the  correlation  functions,  we  make  a  quasinormality 
assumption  which  expresses  the  fourth-order  moment  appearing  in  Eq. 

(30)  according  to  the  scheme 

<ABCD>  =  <AB>  <CD>  +  <AC>  <BD>  +  <AD>  <BC>,  (32) 

valid  if  the  four  field  components  would  approximate  a  joint  normal 
distribution.  This  approximation  does  not  alter  the  homogeneity  of  the 
system,  which  ultimately  causes  the  1/f  spectrum.  This  approximation 
yields  the  closed  equation 

dwap(k,x)/dx  +  pk2wap(k,x) 

=(2n/L)3waj(k,T)Sk,wirT1(k',0)Rjmip(l<,k').  (33) 

Isotropic  turbulence  requires  wap  =  Ai(k)5ap  +  A2(k)kakp,  with  coefficients 
Ai  and  A2  related  through  Eq.  (26),  yielding 


wap(k,t)=  (1/2)[5ap  -  kakp/k2]u(k,x) , 


(34) 


where  u(k,x)  =  ZaWaa(k,t).  Therefore,  the  scalar  correlation  function  u(k,x) 
satisfies  the  dynamical  equation  of  homogeneous,  isotropic,  stationary 
turbulence 


3v(k,x) 

d|x| 


+  v(k,x)  =  -  £v(k,x)J^f  r|y[1  -(l^V)2]v(k',0)I 


(35) 


where  v(k,x)=k-3u(k,x),  and  x  =  pxk2  is  a  dimensionless  variable  replacing  t. 
We  convince  ourselves  that  the  integral  is  independent  of  k,  provided 
v(k,x)  does  not  depend  on  its  first  argument,  by  setting  k'/ksK.  This  yields 
a  solution.  However,  with  v=e’mixl  we  get  a  logarithmic  divergence  at 
k=0.  We  look  for  an  exact  solution  of  the  form4-7 


v(k,x)  =  hkee'lxlm(k)>  0r  u(k.x)  =  (h/k3-e)  e-|x|m(k)f  (36) 

where  m(k)  is  very  close  to  a  constant,  almost  independent  of  k,  and  h  is  a 
constant  proportional  to  the  intensity  of  the  turbulence,  or  the  turbulence 
level.  Substituting  this  into  Eq.  (35)  and  performing  the  integration,  we 
obtain  a  finite  result  only  for  0<e<2: 

m(k)  =  1  +  h  r(e)ke,  with  r(e)  =  [2rc2cotan(E7i/2)]/[(1  -e2)(3-e2)].  (37) 

We  notice  that  m(k)  is  indeed  practically  constant  when  0<e«1  is  very 
small,  arbitrarily  small.  The  value  e=0  leads  to  a  logarithmic  divergence, 
but  we  can  set  e=0  for  practical  purposes. 

The  spectral  density  corresponding  to  Eq.  (36)  with  e=0  is 


wap(co)  =  (1  /rc)  Jcoscoxdx  Jwap(k,x)d3k 
0 


CO  CO 

=  3  5ap  Jk2dk  Jdx  u(k,x)coscox  =  “  h5ap 
0  0 


*  mkdk  _  E  h 
co2+m2k2  _  3  o)  “P- 


(38) 
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This  is  a  1/f  spectrum.  At  the  low  frequency  end  we  do  not  get  a 
divergent  spectral  integral,  because  the  more  exact  form  of  the  spectrum 
w:th  a  finite  small  e«1  is4*7 


Vnk1+£dk  1  [x1+edx 

co2+m2k4  -  co1-e/2mE/2  1+x4 

0  0 


(39) 


which  is  proportional  to  fe/21.  It  is  interesting  to  note  that  for  e<<1 
cotan  £7t/2  -  2/erc  in  Eq.  (37),  and  that  the  value  of  e  calculated  from  Eq. 

(37)  is  therefore  proportional  to  h,  or  to  the  intensity  of  the  turbulence. 
This  feature  of  the  classical  theory4'7  is  expressed  with  fascinating 
clarity  in  the  quantum  form  of  the  theory,  where  £  is  replaced  by  2aA 
which  also  appears  as  a  intensity  factor  multiplying  the  quantum  1/f 
noise. 

The  essential  element  which  led  to  the  1/f  spectrum  in  the  classical 
turbulence  theory  is  the  nonlinearity  of  the  equations  of  motion,  caused  by 
the  reaction  of  the  electric  currents  back  on  themselves  via  the  generated 
electromagnetic  field.  The  same  feedback  reaction,  via  the 
electromagnetic  field,  also  caused  the  nonlinearity  in  the  quantum  1/f 
theory,  and  in  QED  in  general,  leading  in  the  same  way  to  an  identical  1/f 
spectrum,  this  time  with  a  physically  more  meaningful  e= 2aA.  This 
nonlinearity  induces  the  coupling  between  various  scales  of  turbulence 
and  leads  to  the  dynamical  equilibrium  between  eddies  of  all  sizes, 
expressed  by  the  1/f  spectrum.  In  the  e=0,  or  aA=0,  limit,  this  dynamical 
equilibrium  assumes  both  for  the  quantum  case  and  for  homogeneous, 
isotropic,  turbulence  in  the  unbounded  semiconductor  sample  the  simplest 
form,  characterized  by  scale-homogeneity,  or  scale  invariance.  Indeed, 
replacing  for  e=0  in  Eqs.  (35)  and  (36)  k  and  k’  by  Xk  and  Xk\  while  also 
replacing  x  by  x/X2,  (or  o  by  X2co),  Eq.  (35)  is  not  affected,  and  X  drops  out. 
We  conclude  that  in  the  weak  turbulence  limit  (e=0)  we  obtain  perfect 
self-similaritv  of  the  turbulence  process  at  all  scales  in  space  and  time, 
classically  and  quantum-mechanically.  The  implied  scale  invariance  is 
caused  by  the  absence  of  any  characteristic  length  or  time  scale,  or  by  the 
presence  of  a  sliding  scale.  Indeed,  the  frequency  scale  (ik2  is  a  function 


of  the  size  of  the  eddies,  given  by  the  wave  number  k  which  can  have  any 
value.  The  actual  frequency  and  wave-number  spectra  are  closely  related 
fractals,  but  in  the  weak-turbulence  limit  they  approach  an  exact  1/f  and 
1/k3  spectrum  respectively.  In  fact,  we  are  here  understanding  the 
nonlinear  dynamics  which  shapes  this  fractal  for  the  first  time. 

III. 2  Turbulence  Theory  for  Drude  Electrons  in  Metals 
Our  classical  turbulence  theory  can  be  extended  to  the  case  of 


metals  or  degenerate  extrinsic  semiconductors7  in 
system  of  integro-differential  equations  is  quite 

the  Drude  model.  The 
different, 

w+  =  -eE  -  (e/cfvxB  -  (1/n)VP 

(40) 

VxE  =  ~(1/c)3B/at 

(41) 

VxB  =  -(47cen/c)v 

(42) 

V-B  =  0, 

(43) 

and  leads  to  a  third-order  nonlinearity7  in  the  resulting  closed  equation  of 
motion,  or  nonlinear  field-equation,  which  replaces  Eq.  (28): 

3B(k,t)/3t  +  vk2B(k,t)  =  -(cMrcneJkxJdSk'Btk-k’.tMk'xBfk’.t)].  (44) 

This  is  again  in  the  form  of  Eq.  (3),  with  p=2  in  Eq.  (4).  We  thus  expecta 
1/f  spectrum  in  this  system  as  well.  This  time  we  only  sketch  the 
derivation.  The  corresponding  infinite  chain  of  equations  for  the 
correlation  tensors  now  goes  in  steps  of  one.  As  was  shown  above,  it 
went  in  steps  of  two  for  semiconductors.  The  third-order  correlation  can 
be  eliminated  between  the  first  and  second  equations  in  the  chain.  The 
resulting  dynamical  equation7  for  homogeneous,  isotropic,  stationary 
turbulence,  which  replaces  Eq.  (35),  with  the  same  notations,  using  e3  as 
the  unit  vector  of  the  third  axis,  is 


(45) 


92v(k,x) 
9x2  + 


v(k,x) 


1  rd3K  1  +<'63 

2  k3  |e3+K|3 


(j£^-]v(Kk,K2x)v(k|e3+K|)x|e3+K|2). 


This  also  admits,  in  the  e=0  limit  of  weak  turbulence,  a  solution 
v(k,x)  which  does  not  depend  on  the  first  argument,  and  u(k,x)  =  k'3e- 
xm(x),  this  time  with  an  x-dependent  m.  With  the  change  of  variables  t=coT 
and  k'=k/Vco  in  the  second  (middle,  involving  u)  form  of  Eq.  (38),  x  remains 
invariant,  and  a  factor  1/cd  will  appear  in  front  of  the  integrals  which 
themselves  will  just  yield  a  constant  factor  independent  of  co,  x  or  k.  W  e 
thus  obtain  again  a  universal  1/f  spectrum.  As  is  shown  in  detail 
elsewhere7,  this  1/f  spectrum  is  expressed  in  the  corresponding  current 
and  voltage  fluctuations  which  can  be  observed  in  the  semiconductor  or 
metallic  medium.  We  conclude  that  the  1/f  spectrum  is  a  general  property 
of  electrically  conducting  systems  in  interaction  with  the 
electromagnetic  field,  a  property  which  is  caused  by  the  nonlinearity  of 
the  system  of  carriers  and  field  in  mutual  interaction  due  to  the  absence 
of  a  characteristic  scale  in  the  nonlinear  equation  of  motion,  and  which 
finds  its  clearest  expression  in  the  Quantum  1/f  Effect. 


III. 3  Theory  of  Highway  Traffic  Fluctuations 

Musha  and  Higuchi8-9  discovered  the  1/f  spectrum  of  highway  traffic 
fluctuations  empirically  in  1977,  also  developing  a  model  based  on  a 
postulated  linear  dependence  of  the  average  traffic  speed  v  on  the  linear 
concentration  of  cars  n(x)  on  the  road,  v=v0(1-n/ns).  The  model  was 
recognized  to  be  similar  to  Burger's  model  of  turbulence,  and  was 
simulated  numerically  leading  to  a  1/f-like  spectrum  at  low  frequencies. 

The  present  paper  develops  a  statistical  turbulence  theory  for  the 
Musha  model,  showing  how  in  the  low  wave  number  and  low  frequency 
limit  a  universal  1/f  spectrum  emerges.  In  this  limit,  this  is  a  particular 
case  of  the  author's  sliding-scale-invariant  ciass  of  nonlinear  systems, 
all  characterized  by  the  presence  of  a  universal  1/f  spectrum.  The 
author's  classical10-11  and  quantum12  1/f  theory  is  another  example  in 
the  same  class13. 
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Musha  writes  the  traffic  current  J  as  a  sum  of  drift  and  diffusion 
currents 

J  =  nv  -  D3n/3x  =  nv0  -  (v0/ns)n2  -  D3n/3x,  (46) 

where  D  is  a  "diffusion  coefficient".  The  equation  of  continuity  is8-9 

0  =  3n/3t  +  3J/3x  =  (3/3t  +  v03/3x)n  -  2(v0/ns)n3n/3x  -D32n/3x2.  (47) 

This  is  Musha’s  fundamental  equation  of  traffic  dynamics,  which  also  was 
written8  in  a  system  of  reference  defined  by 

X'  =  -X  +  V0t,  r  =  t  (48) 

in  the  final  form8 

3n/at*  +  an3n/3x'  =  D32n/3x'2,  (49) 

where  a  =  2v0/ns. 

We  normalize  the  concentration  to  ns/2,  the  coordinate  along  the 
road  x  to  D/v0,  and  the  time  to  D/v02,  thereby  obtaining  the  dimensionless 
form  of  the  fundamental  traffic-dynamical  equation 

3n(x,t)/3t  +  n(x,t)3n(x,t)/3x  =  32n(x,t)/3x2,  (50) 

where  we  did  not  bother  to  change  the  notation,  and  returned  to  the 
original  unprimed  notation.  Expanding  in  a  Fourier  series  over  the  interval 
L,  we  obtain 

3n/3t  +  iLk-k'n(k')n(k-k')  =  -k2n(k,t).  (51) 

Comparing  Eq.  (50)  with  Eqs.  (1),(2),  or  Eq.  (51)  with  Eqs.  (3), (4  ,  we 
notice  that  this  time  our  criterion  is  not  satisfied,  due  to  the  r.h.s.  term. 
However,  in  the  limit  of  small  k  that  term  becomes  negligible,  and  thus, 
our  criterion  becomes  applicable,  and  we  should  get  a  1/f  spectrum  in  the 


low  frequency  limit.  This  would  correspond  to  p=1  in  Eqs.  (2)  and  (3).  We 
derive  this  1/f  spectrum  below. 

Defining  the  autocorrelation  function 

A(^,x)  =<n*(x,t)n(x+^,t+x)>  =Ek<n*(M)n(k,t+x)>eikS  =LkU(k,x)elk^,  (52) 

as  a  turbulent  ensemble  average  <...>,  in  quasistationary  homogeneous 
conditions  we  obtain  by  multiplication  of  Eq.  (51)  with  n*(k,  t-x) ,  after 
ensemble  averaging, 

3U/3x  +k2U(k,x)  +  iZk'k'<n*(k,  t-t)n(k,,t)n(k-k',t)>  =  0.  (53) 

Linking  the  complex  function  U(k,x)  to  a  third-order  correlation  function, 
this  is  the  first  equation  of  an  infinite  chain  of  equations  connecting  the 
Nth  order  to  the  order  N+1  correlation  function.  Applying  the  operator 
-3/3x  +k2  to  Eq.  (53),  and  using  the  complex  conjugate  of  Eq.  (51)  to  define 
the  action  of  this  operator  on  the  first  factor  inside  the  averaging 
brackets  in  Eq.  (53),  we  obtain 

-3211/3x2  +  k4U  =  -iEk,k,<[iZk"k"n*(k",t-x)n*(k-k",t-x)]n(k',t)n(k-k,,t)> 

-  Zk',k-k'k"[U(k",0)U(k,,0)5k,o  +U(k',i)U(k-k\t)8k-k'  +U(k,,t)U(k-k',t)8k-,k.ie]. 

(54) 

The  last  form  was  obtained  by  approximating  the  fourth-order  correlation 
function  with  its  expression  in  terms  of  the  second-order  correlation 
applicable  for  Gaussian  processes: 

<ABCD>  =  <ABxCD>  +  <ACxBD>  +  <ADxBC>.  (55) 

This  approximation  was  used  by  Heisenberg14  in  his  turbulence  theory,  and 
was  successfully  verified  in  its  practical  applicability  by  Uberoi15  and 
Batchelor16. 

In  the  limit  L-^,  setting  U(k,x)L/27i=  u(kx),  we  write  for  k*0  Eq.  (54) 
in  the  form 


-32u/3x2  +  k4u  =  J[k'2u(k',x)u(k-k',  x)  +  k'(k-k')u(k,,x)u(k-k',x)]dk’ 


or 

-32u(k,t)/3x2  +  k4u(k,x)  =  kjk'uCk'.tMk-k',  t)dk\  (56) 

This  is  our  fundamental  dynamical  equation  of  traffic  turbulence.  All 
unspecified  limits  on  integrals  are  from  -<»  to  as  we  mentioned  earlier. 
From  the  symmetry  A(^,x)  =A(-£,-t)  =A*(4,t)  we  see  that  a  physically 
acceptable  solution  of  this  equation  must  satisfy  the  conditions 

u(k,t)  =  u(-k,-t)  =  u*(-k,-t)  =  u*(k,-x).  (57) 

A  solution  of  the  form 

u(k,t)  =  V(k)e"imkT,  with  V(k)  =  V(-k)  =  V*(k)  (58) 

where  m  is  a  real  constant,  substituted  into  Eq.  (56),  yields 

k(m2  +  k2)V(k)  =  Jk'V(k,)V(k-k')dk'.  (59) 

In  the  low  wave  number  region  k«1  (i.e.,  k«v0/D)  we  neglect  the  k2  term 
which  arises  from  diffusion,  and  get 

m2kV(k)  =  Jk'V(k,)V(k-k,)dk’.  (60) 

A  solution  V(k)  =  Clkl6-1 ,  with  arbitrarily  small  e>0,  yields  the  value  of  C, 
independent  of  k  only  in  the  limit  of  e-»0: 


m2/C  =  Jk'lk’l^-i |i -k7k|e-idk7k  =|k|E  Jk|k|e-1|1  -KpidK 

-oo  -oo 

2 

=  |k|e  ||K|e|1-K|£-1o,K  «  2/e  =  I  «  const, 
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(61) 


where  we  used  the  substitution  kVk=K.  The  divergence  at  e=0  disappears 
when  we  return  to  a  finite  value  of  the  maximal  road  length  L,  which 
corresponds  to  a  minimal  k  (or  k')  value  k0=27i/L,  and  transforms  the 
integrals  back  into  sums. 

Eq.  (61)  establishes  a  proportionality  between  the  level  of  the 
turbulence, described  by  C,  and  the  magnitude  of  the  small  parameter  e  as 
in  earlier  turbulence  calculations10'11.  The  same  fundamental  feature  is 
more  clearly  expressed,  without  the  pseudo-singularities  present  here,  in 
the  quantum  1/f  theory12,  where  e  becomes  the  infrared  exponent  aA 
known  from  quantum  electrodynamics.  Indeed,  aA  is  present  there,  just  as 
we  see  it  come  in  here,  both  as  a  factor  in  front  of  the  final  result,  and  as 
a  small  defect  in  the  exponent  of  the  frequency.  We  also  mention  that  the 
apparent  absolute  determination  of  C  by  Eq.  (61)  is  just  an  artifact  which 
reminds  us  that  we  omitted  a  source  term  in  Eqs.  (47),  (49-51),  (53-54), 
(56),  which  comes  in  as  a  8  function  of  time,  h(k)5(t)  on  the  r.h.s.  of  Eq. 

(53),  only  if  we  assume  the  excitation  in  Eq.  (47)  to  depend  on  time  like 
white  noise.  Once  that  source  term  is  written  explicitly,  it  will  allow  us 
to  determine  C  as  this  was  done  earlier11  for  a  different  equation,  and  it 
does  not  affect  our  equations  for  x  *  0. 

We  can  rewrite  Eq.  (52)  in  the  form 

A(£,x)  =  <n‘(x,t)n(x+^,t+T)>  =  EkU(k,x)eikS  =  Jdk  u(k,x)eik^ 

=  cje'k(£-mx)|k|e'i dk  =  |^-mx|'E(m2/l)JeiK|K|E-1dK  *  m2|^-mx|-E.  (62) 

Here  we  have  used  the  substitution  k|£-mx|=K,  and  we  have  taken  the  limit 
e->0  in  the  integral  and  in  I,  noting  that  they  exhibit  the  same  divergence, 
and  therefore  can  be  simplified  in  the  limit.  The  cancellation  of  the 
divergences  in  the  expression  of  the  autocorrelation  function  shows  that 
this  important  function  is  finite  even  in  the  continuum  limit.  It  is  only 
the  spectral  distribution  of  turbulence  which  exhibits  a  singularity  at  low 
wave  numbers  and  frequencies. 

Since  only  the  limit  e=0  satisfies  the  fundamental  dynamical 
equation  of  traffic  turbulence,  Eq.  (62)  means  that  the  autocorrelation 
function  is  a  constant.  It  also  indicates  that  the  constant  has  to  be 
interpreted  as  the  limit  of  a  slowly  decreasing  power  law  which  depends 


only  on  |£-x|;  therefore  the  spectrum  in  wave  numbers  is  the  same  as  the 
spectrum  in  frequencies.  The  Fourier  transform  with  respect  to  x  is  the 
spectral  density 

Sn’^.co)  =|<n*(x,t)n{x+^t+t)>e'^dx  =  c||(dk/|k|i-£)e'k(^-mt)+icoTdT.  (63) 

To  leave  the  moving  (primed)  frame  and  return  to  the  system  of  reference 
at  rest,  we  replace  co  by  co-kv0  in  the  last  integrand,  actually  just  by  co-k  in 
our  (v0=1)  normalization: 

Sn(£,co)  =j<n*(x,t)n(x+^,t+x)>e'(w_kv,,)'cdx 

=cJJ(dk/k1e)eik(^-m  x)+i(co-k)xc|i;.  (64) 

The  spectral  density  of  concentration  fluctuations  in  a  given  point  is 
obtained  by  setting  £=0.  The  integration  with  respect  to  x  yields  a  delta 
function  27c8[(o-(m+1)k],  and  we  finally  obtain  the  spectrum 

Sn(0,co)  -  C[(m+1)/co]l-e  (65) 

where  again  the  constant  will  turn  out  to  be  finite  when  we  coarse  grain 
the  integrals.  As  mentioned  above,  the  constant  will  actually  be 
determined  by  the  excitation  term  omitted  in  Eq.  (56). 

According  to  Eq.  (2)  we  can  write  J(k,co)  =  (co/k)n(k,co)  in  the  system 
at  rest.  Therefore,  for  J  we  include  a  factor  (co/k)2  into  the  integrand  of 
Eq.  (64),  and  we  get  again  the  spectrum 

Sj(^,co)=/<J*(x,t)J(x+4,t+x)>e'(w-kv.)X(jx 

=  cJJ(co2dk/k3'e)eik(^— mx)+i(co-k)TCjx;  (66) 

Sj(0,co)  =  C(m+1)2  [(m+IJ/coj1*6  =  Const/co.  (67) 

Due  to  the  neglect  of  the  k2  term  in  Eq.  (59),  the  universal  1/f 
spectrum  will  be  limited  towards  high  frequencies  by  cd<v02/D,  and  the 
1/k-spectrum  by  k<v0/D.  This  limitation  occurs  because  only  Eq.  (60) 
satisfies  the  condition  of  being  free  of  any  characteristic  scale,  thus 
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exhibiting  our  sliding  scale  invariance13,  while  Eq.  (59)  does  not  have  this 
property.  We  have  thus  constructed  a  statistical  dynamic  theory  of 
traffic  turbulence,  proving  analytically  Musha's  earlier  result,  without  any 
pretension  of  mathematical  rigor.  Traffic  turbulence  arising  from 
instabilities  of  the  laminar  traffic  flow  can  be  considered  as  a  form  of 
classical  fluid-dynamical  chaos. 

11 1.4  Quantum  1/f  Noise  (QED) 

The  nonlinearity  causing  the  1/f  spectrum  of  turbulence  in  both 
semiconductors  and  metals  is  the  reaction  of  the  field  generated  by 
charged  particles  and  their  currents  back  on  themselves.  The  same 
nonlinearity  is  present  in  quantum  electrodynamics  (QED),  causes  the 
infrared  divergence,  infrared  radiative  corrections  for  cross  sections  and 
process  rates,  and  causes  the  quantum  1/f  effect.  The  basic  equations  of 
QED  are 

a^3^Av=Jvi  Jv  =  -je\jfyv\jr,  (68) 

(3^  +  m)Y  =  -ieAvyv\j/  (69) 

Here  3^  is  a  shorthand  for  3/3x^,  yM-  are  the  Dirac  matrices,  the  Greek 
indices  run  from  0  to  3,  and  y  is  the  Dirac  spinor  field. 

We  can  formally  solve  Eq.  (68)  for  the  potential  Av  in  momentum 
space  (i.e.  by  Fourier  transforming) 

Av  =  (1/ana^jv,  (70) 

which  can  be  substituted  into  Eq.  (70),  yielding  with  the  expression  (68) 

of  the  current  Jv, 

(3^  +  mc/h)\)/  =  -tx[(1/3^3|i)\(/yvx|/]yv'V/-  (71) 

The  integral  operator  (1/3^3^)  has  a  kernel  1/q2,  where  q  is  the  Fourier 
four-vector  corresponding  to  x.  Eq.  (71)  contains  a  third-order 
nonlinearity  and  is  therefore  similar  to  the  dynamical  equation  (29) 
describing  turbulence  in  semiconductors.  In  the  spinless  non-relativistic 
limit,  Schrodinger's  equation  yields 
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ay/at  =  [-ihV  +  A]2y 


(72) 


or,  with  Eq.  (68), 

2mihay/at  =  [-ihV  -  i(e2h/mc2)(1/3na4)y*Vy]2y  (73) 

At  very  low  frequencies  the  last  term  in  rectangular  brackets  is  dominant 
on  the  r.h.s.,  leading  to 

2iay/at  =  -h[(92/mc2)(1/aua^)y*Vy]2y  (74) 

Both  Eq.  (71),  in  which  we  neglect  at  low  wave-numbers  the  terms  in 
a/axj  (j=1,2,3)  and  eliminate  the  term  with  the  rest  mass  m  through  a 
substitution  \j/=<t>eimTt>  and  right  away  Eq.  (74),  satisfy  our  criterion  with 
p=2  in  Eq.  (2).  Therefore,  we  expect  a  1/f  spectrum  of  current 
fluctuations,  i.e.,  of  cross  sections  and  process  rates  in  physics.  This  is 
in  agreement  with  the  well-known,  and  experimentally  verified,  results  of 
the  conventional  Quantum  1/f  Theory. 

In  conclusion,  we  realize  that,  both  in  classical  and  quantum 
mechanical  nonlinear  systems,  the  limiting  behaviour  at  low  wave 
numbers  is  usually  expressed  by  homogeneous  functional  dependences,  (as 
shown  in  Eq.  2),  leading  to  fundamental  1/f  spectra  on  the  basis  of  our 
criterion. 


IV.  ANALYTICAL  CALCULATION  OF  MOBILITY  FLUCTUATIONS 
IN  SEMICONDUCTORS,  BASED  ON  THE  QUANTUM  1/f 
CROSS-CORRELATION  FORMULA 

Together  with  the  graduate  student  Thomas  H.  Chung,  the  author  has 
performed  an  analytical  calculation  of  mobility  fluctuations  in  silicon  and 
gallium  arsenide,  using  the  new  quantum  1/f  cross-correlations  formula 
derived  by  the  author  in  the  previous  AFOSR  Grant  period,  and  included  in 
the  July  1989  Final  Technical  Report.  This  calculation  is  of  mjor 
importance  for  the  1/f  noise-related  optimization  both  of  the  two  types 
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of  materials,  and  of  the  many  devices  constructed  with  them  for  military 
and  civilian  applications  in  the  electronic  and  opto-electronic  industry. 

The  new  cross-correlation  formula  gives  the  cross-spectral  density 
which  describes  the  way  in  which  symultaneous  quantum  1/f  scattering 
rate  fluctuations  AW  obseved  in  the  direction  of  the  outgoing  scattered 
wave-vector  K'  are  correlated  with  those  in  the  K"  direction,  when  the 
two  corresponding  incoming  current  carriers  have  the  wave  vectors  Ki  and 
K2: 

Saw(Ki  ,K';K2,KH;f) 

=  (2a/37cf)(h/m*c)2WKi,K'W  K2,Km[(K,-Ki)2+(K"-K2)2]5Ki,K2-  (75) 

The  form  conjectured  by  us  earlier  had  2(K'-Ki)(K"-K2)  in  place  of  the 
rectangular  bracket. 


IV. 1  Impurity  Scattering 

For  impurity  scattering  of  electrons  in  solids,  fluctuations  At  of  the 
collision  times  t  will  cause  mobility  fluctuations 

Apband(t)  -  [e/m*«v2»]ZKVK2At(t)nK,  (76) 

where  «v2»  is  both  the  average  over  all  states  of  wave-vectors  K,  with 
occupation  numbers  nK,  in  the  conduction  band,  and  the  thermal 
equilibrium  average  of  the  quadratic  carrier  velocities.  With  the  help  of 
the  relation 

1/t(K)  =  (V/87t3)J(l  -  cos07cos0)WK,K'd3K',  (77) 

the  mobility  fluctuations  are  reduced  to  fluctuations  of  the  elementary 
scattering  rates  Wk,k’.  governed  by  Eq.  (75).  Here  V  is  the  volume  of  the 
normalization  box  which  disappears  in  the  final  result,  0  and  0' 
respectively  the  angles  K  and  K'  form  with  the  direction  of  the  applied 
field.  One  finally  obtains  after  tedious  multiple  integrations 

F'2SAn(f)  =  [2567taK2e4hi2/3m*8Z4e8Ni2](1/f)IKK10[ln(1+a2)- 
a2/(1  +a2)]-3[(2a2+a4)/(1  +a2)-2ln(1  +a2)]F(EK)[EKVK2t<K)F(EK)]-2,  (78) 
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where  a=2K/»c,  K2=e2n(T)/ekBT,  n(T)  is  the  electron  concentration, 
F(EK)=exp(Ep-EK)  for  non-degenerate  semiconductors,  Nj  the  concentration 
of  impurities  of  charge  Ze  and  e  the  dielectric  constant.  The 
corresponding  partial  Hooge  parameter  for  impurity  scattering  is  thus 


ocj  =  [4V2mKh5Nc/3m*7/2(kBT)3/2c2]  JdxxH/2e-x 

0 

[ln(bx+1  )-bx/(bx+1  )]-3[(2bx+b2x2)/(bx+1  )-2ln(bx  +  1 )] 

oo 

{  Jdxx3e-X[in(bx+1  )-bx/(bx  +  1 )]-  'J-2-  (79) 

0 


IV. 2  Acoustic  Electron-Phonon  Scattering 
In  this  case  the  calculation  is  similar,  and  leads  to  the  result 


ocac  =  [32KaNcm*C7h3/3c2kBT)4]{(1/R2)  Jdxx~4 

1 

[(x-1)7/7+(R  +  1)(x-1)6/6+R(x-1)5/5] 

[(x-1  )5/5+(R+1 )  (x- 1  )4/4  +  R(x-1  )3/3]exp(-x2/4R) 

1 

+  Jdxx'4[(x+1  )5/5-(x+1  )6/6+(x-1  )5/5  +  (x-1  )6/6] 

0 

[(x  +  1  )3/3  +  (x-1  )4/4  +  (x-1  )3/3-(x+1  )4/4]exp(-x2/4R) 

oo 

+  Jdxx-4[(x+1)5/5-(x+1)6/6][(x+1)3/3-(x+1)4/4]exp(-x2/4R)},  (80) 
1 

where  R=kBT/2m*Ci2,  Ci  is  the  deformation  potential,  and  Nc  is  the 
effective  density  of  states  for  the  conduction  band. 
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IV. 3  Non-Polar  Optical  Phonon  Scattering 


This  time  one  obtains 

oo 

an.o.ph  =  [8jt\/ 2hoJoaNch2/3m‘5/2c2co0]{  Jdxx5/2 

0 

[(F+1)(X-1)1/20(x-1)  +  F(x+1)1/2]-4 

[(F+1  )2(x-1  )(2x-1)0(x-1  )+F2(x+1  )(2x+1  )]exp(-ha)0x/kBT)} 

oo 

|dxx3/2[(F+1  )(x-1  )1/2e(x-1  )  +  F(x+1  )1/2]-iexp(-hcooX/kBT)}'2,  (81) 


where  F=[exp(hco0/kBT)-1  ]_1 ,  and  w0  is  the  optical  phonon  frequency. 


V.  CALCULATION  OF  COLLECTOR  1/f  NOISE  IN  BJTs 

Based  on  a  consequent  application  of  the  conventional  quantum  1/f 
approach,  we  show  that  for  the  calculation  of  conventional  quantum  1/f 
noise  in  the  collector  current  of  BJTs  the  effective  lifetime  x  of  minority 
carriers  in  the  base  must  be  used  instead  of  their  diffusion  time  xp 
through  the  base.  This  allows  us  to  represent  the  effective  number  of 
carriers  in  the  denominator  of  our  quantum  1/f  formula  as  Ix/e  for  noise 
in  the  collector  current,  and  the  theory  of  Van  der  Ziel  and  Kleinpenning 
remains  valid  with  xp  replaced  by  x.  A  long-standing  difficulty  present  in 
the  application  of  bulk  noise  sources  to  the  collector  current  is  thus 
removed  in  a  natural  way,  and  the  agreement  of  the  quantum  1/f  theory 
with  the  experiment  is  substantially  improved,  particularly  for  narrow 
base  BJTs. 

1/f  noise  in  bipolar  junction  transistors  (BJTs)  was  elegantly 
treated  by  van  der  Ziel17'19  who  applied  a  Hooge-type  approach  similar  to 
Kleinpenning’s  treatment20  of  pn  junctions,  and  used  experimental  data  to 
determine  the  Hooge  constant  which  was  in  turn  compared  with  the 
quantum  1/f  theory.  However,  since  the  BJT  is  a  minority  carrier  device, 
it  requires  the  application  of  the  quantum  1/f  (Handel)  equation1 9. 1 -3, 21 
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from  the  beginning,  for  the  correct  interpretation  of  the  number  of 
carriers  in  the  denominator  of  the  Langevin  noise  source. 

In  the  most  elementary  model22  of  a  BJT,  the  collector  current  lc 
arises  from  minority  carriers  injected  from  the  emitter  into  the  base, 
which  diffuse  across  the  width  Xb  base  and  are  then  all  swept  across  the 
reverse-biased  collector  junction  by  the  built-in  field  of  the  junction.  If 
we  neglect  the  usually  small  leakage  current  of  the  collector  junction  and 
the  small  fraction  of  the  carriers  recombining  in  the  base,  we  get  for  a 
n+pn  BJT 

lc  =  AqDn[n0Bexp(qVBE/kT)/XB],  (82) 

where  A  is  the  cross  sectional  area  of  the  base,  q=-e  is  the  charge  of  the 
minority  carriers  in  the  base,  Dn  their  diffusion  coefficient  in  the  base, 
nB(0)  =  nOBexp(qVBE/kT)  is  the  electron  concentration  at  the  limit  of  the 
emitter  space  charge  region,  Vbe  is  the  applied  base  -  emitter  voltage,  and 
Xb  is  the  width  of  the  base.  The  expression  in  rectangular  brackets  is  the 
electron  concentration  gradient  calculated  with  the  boundary  condition  of 
a  vanishing  electron  concentration  at  the  limit  of  the  collector  space 
charge  region.  We  assume  the  base  to  be  much  narrower  than  the  electron 
diffusion  length  Ln  =  (Dn/t),  XB<<Ln,  but  sufficiently  wide  to  avoid  ballistic 
electron  transport  across  the  base.  Usually  Xb  is  a  fraction  of  a  micron. 

Quantum  1/f  fluctuations  of  the  collisional  cross  sections  of  the 
electrons  in  the  base  will  yield  fluctuations  of  the  diffusion  constant,  and 
of  the  mobility  (5Dn/Dn=8ji/p) 

51c  =  Aq(5Dn)[n0Bexp(qVBE/kT)/XB].  (83) 

The  corresponding  spectral  density  of  fractional  fluctuations  I'2S|C  is 

IC'2<(8lc)2>f  =  Dn‘2<(5Dn)2>  =  p‘2<(5p)2>  =  ccn/fN.  (84) 

In  the  last  step  our  quantum  1/f  equation19-1 -3  was  used,  where  N  is  the 
number  of  carriers  which  define  the  scattered,  or  diffused,  current 
leaving  the  base  and  emerging  in  the  collector,  while  an  =aAn  is  the 
effective  quantum  1/f  noise  coefficient,  or  Hooge  constant.  The  number  of 
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electrons  N  is  thus  determined  by  the  effective  lifetime  x  of  the 
electrons,  which  will  be  slightly  lower  than  the  lifetime  in  the  unbounded 
collector  material,  due  to  the  collector  lead  contact  processes,  and  due  to 
lateral  surface  recombination.  Indeed,  we  can  write  N  =  tlc/q.  Thus  we 
finally  obtain  the  spectral  density  of  the  collector  current  fluctuations 

S|C  =  cxnlq/fx,  (85) 


in  which  x  is  the  effective  lifetime  of  the  majority  carriers  in  the 
collector.  This  expression  is  simpler,  but  similar  to  the  expression 
derived  earlier,  with  the  important  difference  that  now  we  have  a 
lifetime  of  the  carriers  in  the  denominator,  while  before  it  was  the 
usually  much  smaller  diffusion  time  tD  =  Xe2/Dn  of  the  electrons  in  the 
base.  Eq.  (85)  also  implies  that  in  narrow-base  BJTs  of  various  base- 
widths  an  will  be  constant,  as  in  other  devices,  rather  than  an/xD- 
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